Grain-boundary grooving and agglomeration of alloy thin films 
with a slow-diffusing species 



Mathieu Bouvillcl and Dongzhi Chi 

Institute of Materials Research and Engineering, Singapore 117602 

David J. Srolovitz 

Department of Physics, Yeshiva University, 
New York, NY 10033, USA 
(Dated: February 6, 2008) 

We present a general phase-field model for grain-boundary grooving and agglomeration of poly- 
crystalline alloy thin films. In particular, we study the effects of slow- diffusing species on grooving 
rate. As the groove grows, the slow species becomes concentrated near the groove tip so that further 
grooving is limited by the rate at which it diffuses away from the tip. At early times the dominant 
diffusion path is along the boundary, while at late times it is parallel to the substrate. This change 
in path strongly affects the time-dependence of grain boundary grooving and increases the time to 
agglomeration. The present model provides a tool for agglomeration- resist ant thin film alloy design. 
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A common failure mode in polycrystalline thin films 
is grain-boundary grooving through the thickness of the 
film. The grain-boundary shrinks while the surface or 
substrate/film interface extends in order to reduce the 
total interface energy (see Fig. [1]). At equilibrium, 
the surface and interface shapes must correspond to a 
constant chemical potential. For an isotropic interface 
energy, this means constant curvature, i.e. a circular 
arc in two dimensions (2D). The interface thus evolves 
towards a constant curvature shape with a deepenin 
grain-boundary groove (subject to fixed groove angle 
[2[ and grain volume constraints). This can bring the 
surface in contact with the substrate, leading to film 
agglomeration as shown in Fig. 1(d) Many poly- 
crystalline metal silicide films fail in this manner during 
post-reaction annealing 0, B] . 

Grain-boundary grooving has received a great deal of 
attention over the past half-century [l], 0, H 0, @, EJ 

However, the extant models for grain-boundary 
grooving and thin film agglomeration are too idealized 
to be useful to predict agglomeration in most technologi- 
cally interesting materials, such as multicomponent alloy 
films. In particular, few models consider the complexity 
of real grain-boundary grooving — simultaneous bulk, 
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FIG. 1: (color online) The time evolution of the film morphol- 
ogy for h = 10 nm and r = 1. (a) t = 0.15 s, (b) t = 7.5 s, 
(c) t = 63 s, and (d) t = 64.5 s. 



grain-boundary, interface, and surface diffusion; phase 
transformations; etc. In this Letter, we use the phase- 
field method to study grain-boundary grooving and the 
agglomeration of a polycrystalline alloy thin film where 
each alloy component has a unique diffusivity and a phase 
transformation occurs. For concreteness, we choose a 
polycrystalline NiSi film on a silicon substrate as our 
prototype system. We examine the effect of introducing 
a slow- diffusing species (Pt) into the NiSi thin film. 

In our phase-field calculations, the concentrations 
of silicon, nickel, and platinum are csi, CNi, and cp t 
respectively; all three fields are conserved. <fi indicates 
the local fraction of the silicon phase present and ipi the 
local fraction of grain i of the compound silicide phase; 
these are not conserved. The free energy of the system 
is given by [l3[ 
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In Eq. (pp), the integral is over the entire volume of 
the system, V. The first three terms are double-well 
potentials for ipi and 0, which ensure that in the bulk 
ipi and cj) are either zero or one. The fourth and fifth 
terms ensure that at any point in the system, only a 
single phase is present. The third line of Eq. (pQ) fixes 
the stoichiometry (in the substrate Cm = and Civ = 1, 
whereas in the silicide film Cm = Civ = 1/2), while 
the fourth line describes the configurational entropy 
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FIG. 2: The depth of the groove for h — 10 nm as a function 
of time (log-log scale). •: r = 1, M: r = 3 X 10" 2 , □: r = 
10" 2 , o: r = 3 x 10~ 3 , and +: r = 10~ 3 . The curves are offset 
by a half decade, for clarity. The lines are power law fits to 
Eq. J4} with parameters as shown in Table |U 

associated with the metal sublattice of the silicide. The 
gradient terms in the last two lines enforce smooth 
interfaces. The parameters W 7 , X/, Xqb, A, ft, and set 
the overall thermodynamics of the film (Q is the average 
atomic volume). 

The temporal evolution of the compositions and phases 
is related to the free energy G, respectively through the 
Cahn-Hillard equation [14[ 

dc(r,t) m _ SG 
ot oc(r, t) 

and the Allen-Cahn equation [l5[ 

fl^(r,t) _ SG 

dt 6il>(r,t)' [ } 

Here M and L are atom and interface mobilities, respec- 
tively, and S indicates a functional derivative. 

The two-dimensional system consists of two silicide 
grains on a silicon substrate, discretized on a 512 x 96 
mesh. For simplicity, we consider the symmetric case 
where the phase above the film is also silicon (rather 
than vacuum). This allows us to use periodic boundaries 
along both directions, as well as fewer parameters in the 
free energy (e.g. the surface and interface are identical). 
The evolution equations, Eqs. ([2]) and (j3]), are solved 
in reciprocal space using a semi-implicit method, as 
described in Ref. We use the following values for 
the parameters: W = 5, X\ = 30, Xqb = 18, A = 8, 
k B T/n = 0.3, A Si = 20, A Ni = 30, A Pt = 100, ft = 8, 
L = 4, and M = 1. We choose the length and time scales 
to represent a typical case of a 10 nm-thick NiSi film 
on Si, which is experimentally observed to agglomerate 
in around a minute at T — 600 °C. Film thicknesses 
are thus given in nanometers for convenience, not based 
on the interface thickness. One should note that our 
simulations are relevant to grain-boundary grooving and 
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TABLE I: The grooving parameters do (in nm s~ a ) and a 
from Eq. J4j — as determined by fits in Fig. [2] — , t a , the 
time at which the film agglomerates, and c a , the value of 
c at agglomeration for h = 10 nm. For small values of r, we 
quote early and late time values of do and a corresponding to 
the two regimes shown in Fig. [2] t* is the time at which the 
early and late time fits intersect in Fig. [2] 

agglomeration of alloy thin films in general and that 
ratios of times and distances are more important than 
absolute values. 

This Letter focuses on how the introduction of a slow- 
diffusing species into the film affects grain-boundary 
grooving and agglomeration. We perform a series of 
simulations at different ratios of the Pt to Ni mobilities, 
r = Mpt/MNi < 1 (we keep the Ni mobility fixed), 
spanning three orders of magnitude. We note that the 
silicide is a compound in which the Pt substitutes for Ni 
on the metal sublattice. In the present simulations, the 
average concentration of Pt is 7%. 

Figure [2] shows the groove depth d as a function of 
time t for several values of r. The time-dependence of 
the groove depth is of the form 

d = d t a . (4) 

Mullins [2[ found that a — 1/4 if mass transfer occurs 
by surface diffusion and a = 1/2 if by evaporation and 
condensation. At early times the groove depth is well 
described by Eq. (|4|), with a between 0.24 and 0.30 (see 
Table H]), in reasonable agreement with the theoretical 
exponent for the surface diffusion case. The prefactor 
do increases with r; as expected, the presence of a 
slow-diffusing species in the film retards grain-boundary 
grooving. However, Fig. [2] shows that at late times the 
apparent value of a decreases to a lower, constant value. 
This cross-over time t* is larger for smaller r. A change 
of the value of a has also been found in Refs. H and0. 

Clearly, decreasing the mobility of one of the species in 
the film (Pt here) has a major effect on the time required 
for agglomeration to occur. Table [J shows that decreasing 
r from 1 to r = 10 -2 triples the time required for film 
agglomeration, while decreasing r to 10 -3 leads to a 
twenty-fold increase. Although introduction of a slow- 
diffusing species was expected to slow grain-boundary 
grooving and delay agglomeration, the complex interplay 
between morphology and composition in multicomponent 
systems adds unforeseen richness. An alloy with a slow- 
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FIG. 3: (color online) The Pt composition field at t = 9 s for 
/i = 10 nm and several values of r = M-pt/Mm- (a) r = 1, (b) 
r = 10 _ \ (c) r = 10~ 2 , and (d) r = 10~ 3 . 



diffusing species does not behave like a monocomponent 
system with a lower diffusivity, as often assumed. 

Figure [3] shows that for small r there is a notable 
excess of platinum near the grain-boundary /interface 
intersection; the Pt content close to the groove is roughly 
40% larger than the average Pt concentration in the film. 
Indeed the groove can grow only as fast as the metal 
atoms can diffuse away from the deepening groove. When 
Pt diffuses much more slowly than Ni, the grooving rate is 
limited by the time required for the advancing groove to 
'push the Pt out of the way.' On the other hand, if r = 1 
or r = 10 _1 (see Fig. [3]), the Pt composition is nearly 
homogeneous throughout the film. In this case, the 
alloy behaves similarly to a monocomponent system and 
the time evolution of the groove depth follows Eq. (|4]). 
This is consistent with the observations made above for 
the effect of Pt on the rate of grooving (Fig. [2] and 
Table H]) and with experimental observations that show 
that the introduction of Pt into nickel silicide greatly 
slows grooving and agglomeration [171 ]. 

In order to quantify the magnitude of the platinum 
excess observed close to the grain-boundary groove in 
Fig. 02 we measure the maximum Pt content in the 
system as a function of time t, c. Figure [H shows c as a 
function of time for several initial thicknesses. For thick 
films (h = 20 nm), there are two regimes: c initially 
increases and then decreases for t > 12 s (T^j. In the 
thinner films, c exhibits the same early time behavior as 
in the thick films, but at later times the maximum Pt 
concentration grows rapidly. This film thickness effect 
only occurs when the Pt mobility is sufficiently low. 

Since c exhibits the same maximum value at the same 
time in the thicker films (at t « 12 s), this maximum 
cannot be attributed to the finite thickness of the films. 
At early times, grooving is limited by the slow Pt 
diffusion; the groove can grow only as fast as it can 




FIG. 5: (color online) The platinum composition field, cpt, for 
h = 10 nm and r = 10 -3 at (a) t = 30 s and (b) t = 1.7 min. 



'push' the Pt atoms out of the way. This results in a 
Pt build-up close to the groove. As grooving proceeds, 
the driving force for grooving (i.e. the curvature of the 
interface) decreases, so that at later time Pt diffusion is 
able to smooth the Pt concentration profile faster than 
groove advancement can build it. This is what leads to 
the decrease of c around t = 12 s in Fig. [H 

Figure [H shows that in thin films, the behavior of c is 
much more complex than in the thicker films. Following 
the same rise in c as for the thick films at early times, in 
the thinner films, c decreases, then grows again, before 
finally decreasing. This introduces a second maximum 
in c, c max 2, at later times. The first maximum in c, 
c max i, occurs in all thick films at the same time. The 
second maximum, on the other hand, depends on h (for 
sufficiently thin films, it may occur so early that no c max i 
is observed). 

In order to understand how finite thickness can affect 
c, we examine the Pt composition field at different times. 
Along with Fig. |3(d)] F ig. [5] shows cp t for h = 10 nm and 
r = 10 -3 . Figures 3(d) and |5 (b)| correspond to the times 
at which the two c maxima in Fig.[4]occur, while Fig .|5(a) 
corresponds to the intervening minimum. In Fig. [5(a7[ 
the region of large Pt concentration reaches the bottom 
of the film and in Fig. |5(b)[ the Pt composition is nearly 
uniform through the film thickness (but still depends 
on the distance from the grain-boundary, along x). At 
this point, the regions of large Pt concentration near the 
groove can no longer be relaxed by diffusing in the film 
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FIG. 6: (color online) Left and right c y (in %/nm) for h = 
10 nm and r = 10" 3 at (a) t = 9 s, (b) t = 30 s, and (c) t 
1.7 min. The contributions to <4 and c' y from the interface 
and grain-boundary have been removed for clarity. 
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TABLE II: The times, concentrations, and groove depth at 
the two maxima in c vs. time, for h = 10 nm. 



The present simulations demonstrate that introduc- 
tion of a slow- diffusing species can greatly slow grain- 
boundary grooving and agglomeration. The grooving 
leads to pile-up of the slow- diffusing species and the 
grooving rate depends on how quickly this excess can 
be diffusionally relaxed. In the case of grain-boundary 
grooving in thin films, the film thickness is crucial in 
determining the dominant transport mode for the slow- 
diffusing species. When the film is thick compared with 
the groove depth and the slow- diffusing species pile- up, 
transport is dominated by diffusion parallel to the grain- 
boundary. On the other hand, when the film is thin, 
the transport must occur predominantly along the plane 
of the film. In many cases, there is a transition in the 
dominant transport mode and this gives rise to a change 
in the exponent of the grooving power law exponent. The 
time at which this transition occurs depends on the film 
thickness and the mobility of the slow species. 

We gratefully acknowledge S. Hu and L.-Q. Chen for 
their assistance with developing our phase-field simula- 
tion code. This project was funded by the A-Star Visiting 
Investigator Program. 



thickness direction (y-direction), but now must diffuse 
away through the film in the x-direction. This change in 
transport path is very important. 

In order to quantify the driving forces for Pt diffusion 
along x and y, we measure the Pt concentration gradients 
in the x- and y-directions, i.e. — dcp t /dx and c' y = 
dcpt/dy. Figure [6] shows and c' y at the same times as 
Figs. |3(d)| and [51 It confirms that the driving force for 
diffusion along y drops between the two maxima of cp t , 
whereas the magnitude of the driving force for diffusion 
along x varies less over time. While Pt transport is 
dominated by diffusion parallel to the grain-boundary at 
early time for all film thickness, in the thin films, diffusion 
in the y-direction is very limited (by the film thickness) 
at intermediate and late times. This transition in the 
direction of Pt diffusion corresponds to the second rapid 
rise in c in Fig. |5(a)| for the thinner films. At very late 
times, c decreases as the grooving rate slows and Pt has 
sufficient time to diffuse away. 

The times at which the exponent describing the time 
dependence of the groove depth, a, changes at small r 
( cf. Table H]) correlate with the times at which the second 
maximum in c (c max 2) occurs, as shown in Table [TTJ 
There is also agreement for larger values of r; the value 
of a is constant and there is no clear maximum of c. 
Therefore, the change in the grooving exponent a is 
caused by the change in the dominant Pt transport mode. 
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